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Abstract 

Using three- and four-body decays of D mesons produced in semileptonic 6-hadron 
decays, precision measurements of D meson mass differences are made together 
with a measurement of the D° mass. The measurements are based on a dataset 
corresponding to an integrated luminosity of 1.0 fb _1 collected in pp collisions at 
7 TeV. Using the decay D° ->■ K+K-R-1T+, the D° mass is measured to be 

M(D°) = 1864.75 ± 0.15 (stat) ± 0.11 (syst) MeV/c 2 . 

The mass differences 

M(D + ) - M(D°) = 4.76 ± 0.12 (stat) ± 0.07 (syst) MeV/c 2 , 
M(D+) - M(D+) = 98.68 ± 0.03 (stat) ± 0.04 (syst) MeV/c 2 

are measured using the D° — > K + K~tt + tt~ and Df, — > K + K~ir + modes. 
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1 Introduction 



Mesons are colourless objects composed of a quark- ant iquark pair bound via the strong 
interaction. Measurements of meson masses provide observables that can be compared 
to theoretical predictions. For the case of B mesons, precision measurements have been 
reported in recent years by several experiments [ffj3]. In contrast, few precision D meson 
mass measurements exist. 

For the D° mesorj^the current average of M(D°) = 1864.91 ± 0.17 MeV/c 2 , quoted by 
the Review of Particle Physics j3] , is dominated by the measurements of the CLEO [5] 
and KEDR (6) collaborations. Current knowledge of the masses of the D + and Df 
mesons, and the mass splitting between these states, is more limited. The most precise 
determination of the D + mass is made by the KEDR collaboration [6] resulting in 
M(D + ) = 1869.53±0.49 (stat) ±0.20 (syst) MeV/c 2 . In addition, two measurements of the 
mass splitting between the D + and D° mesons by the MRK2 |7| and LGW g collaborations 
have been reported. These are averaged |I] to give M(D + ) - M(D°) = 4.76 ± 0.28 MeV/c 2 . 
No absolute measurement of the -D+ mass with a precision better than the MeV/c 2 level 
exists and the reported values are not in good agreement [4|. More precise measurements of 
the mass difference relative to the D + meson have been reported by several collaborations 
[9}{L3]. These are averaged gjto give M(D+) - M(D + ) = 98.85 ± 0.25 MeV/c 2 . The 



fit of open charm mass data^ leads to M(D+) = 1968.49 ± 0.32 MeV/c 2 . Though this 
value is significantly more precise than the direct measurement, it would still dominate the 
systematic uncertainty on the measurement of the mass in the JB+ — > J/ipDj~ decay 
mode 14 . 

Recent interest in the D° mass has been driven by the observation of the X(3872) state, 



first measured by the Belle experiment 15 and subsequently confirmed elsewhere 16 -20 



This state, with J PC = 1 ++ 21 , does not fit well into the quark model picture, and exotic 



interpretations have been suggested: for example that it is a tetraquark 22 or a loosely 
bound deuteron-like D*°D° 'molecule' 23 . For the latter interpretation to be valid, the 



mass of the X(3872) state should be less than the sum of the D*° and D° masses. Using 
the fitted value of the D° mass and the measured values for the other quantities quoted in 
Ref. |4], the binding energy (Eg) in this interpretation can be estimated to be 

E B = M(D°D*°) - M(X (3872)) 

= 2M(D°) + AM(D*° - D°) - M(X(3872)) 
= 0.16 ± 0.32 MeV/c 2 . 

Therefore, the issue of whether the X(3872) can be a bound molecular state remains open. 
To clarify the situation, more precise measurements of both the X(3872) and D° masses 
are needed. 

In this paper, a measurement of the D° mass using the D° — > K + K~ K~n + decay 
mode is reported. This mode has a relatively low energy release, Q-value, defined as the 
difference between the mass of the D meson and the sum of the masses of the daughter 



lr The inclusion of charge conjugate states is implied. 



1 



particles. Consequently, systematic uncertainties due to the calibration of the momentum 
scale of the detector are reduced. Other four-body D° decay modes are used to provide a 
cross-check of the result. In addition, precision measurements of the D + — D° and Df — D + 
mass differences are made. For the mass difference measurements the D° — > K + K~tt + tt~ 
mode is used, together with the — > K + K~7r + decay, since these modes have similar 
Q- values. 



2 Detector and dataset 

The analysis uses data, corresponding to an integrated luminosity of 1.0 fb _1 , collected in 
pp collisions at a centre-of-mass energy of \fs = 7 TeV by the LHCb experiment during 
2011. The detector response is studied using a simulation. Proton- proton collisions are 



generated using Pythia 6.4 24] with the configuration described in Ref 25 . Particle 



decays are then simulated by EvtGen |26| in which final state radiation is generated 



using Photos 27 . The interaction of the generated particles with the detector and its 



response are implemented using the Geant4 toolkit |28| with the settings described in 



Ref. 29 



The LHCb detector [30] is a single-arm forward spectrometer covering the pseudo- 
rapidity range 2 < rj < 5. It includes a high precision tracking system consisting of a 
silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon-strip 
detector located upstream of a dipole magnet with a bending power of about 4 Tm, and 
three stations of silicon-strip detectors and straw drift tubes placed downstream. The 
polarity of the dipole magnet is reversed at intervals that correspond to roughly 0.1 fb _1 
of coUected data in order to minimize systematic uncertainties. The combined tracking 
system has momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at 
lOOGeV/c, and impact parameter resolution of 20 |xm for tracks with high transverse 
momentum (pt)- Charged hadrons are identified using two ring-imaging Cherenkov 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and pre-shower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identified by a system composed of alternating 



layers of iron and multiwire proportional chambers. The trigger 31 consists of a hardware 
stage, based on information from the calorimeter and muon systems, followed by a software 
stage that applies a full event reconstruction. Samples of open charm mesons produced 
directly in the primary pp interaction (refered to as 'prompt') and in semileptonic decays 
of 6-hadrons are selected by the trigger. Though the prompt sample is larger in size, cuts 
on the decay time of the D meson are applied at the trigger level to reduce the output 
rate. As the reconstructed mass and decay time are correlated, these cuts bias the mass 
measurement. In contrast, no cuts on the D decay time are applied at the trigger level for 
the semileptonic sample, which is therefore used for this analysis. 

The measurements require the momenta of the final state particles to be determined 
accurately. The procedure used to calibrate the momentum scale of the tracking system 



for this study is discussed in detail in Ref. 32 . It is based upon large calibration samples 
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of B + — > J/ipK + and J/ip — > fi + p~ decays collected concurrently with the dataset used 
for this analysis. The use of the large J/ip dataset allows to correct for variations of the 
momentum scale at the level of 10 -4 or less that occur over the course of the data-taking 
period whilst the use of the B + — > J/ipK + allows the momentum scale to be determined 
as a function of the K + kinematics. The accuracy of the procedure has been checked 
using other fully reconstructed B decays together with two-body T(nS) and K® decays. 
In each case the deviation of the measured mass from the expected value is converted 
to an estimate of the bias on the momentum scale (a) taking into account relativistic 
kinematics and QED radiative corrections. The largest value of \a\ found in these studies 
is 0.03 % for the Kg — > 7r + 7r~ decay mode. Conservatively, this is taken as the uncertainty 
on the calibrated momentum scale. This leads to the largest contribution to the systematic 
uncertainty on the mass measurements. 

3 Selection 

The selection uses only well reconstructed charged particles that traverse the entire tracking 
system. All charged particles are required to be within the angular acceptance of the 
spectrometer. This corresponds to 300 mrad in the bending plane of the dipole magnet and 
250 mrad in the orthogonal plane. In addition, the final state particles are required to have 
Pt greater than 300 MeV/ c. Further background suppression is achieved by exploiting the 
fact that the products of heavy flavour decays have a large distance of closest approach 
('impact parameter') with respect to the pp interaction vertex in which they were produced. 
The impact parameter x 2 with respect to any primary vertex is required to be larger than 
nine. Fake tracks created by the reconstruction are suppressed by cutting on the output 
of a neural network trained to discriminate between these and real particles. This cut 
also removes candidates where one of the charged hadrons has decayed in flight. To select 
well-identified kaons (pions) the difference in the logarithms of the global likelihood of the 
kaon (pion) hypothesis relative to the pion (kaon) hypothesis provided by the ring-imaging 
Cherenkov detectors is required to be greater than five (zero). 

Charged particles selected in this way are combined to form D° — > K + K~tt + tt~ , 
D° — > K + K~ K~ir + and Df, — > K + K~ir + candidates. To eliminate kinematic reflections 
due to misidentified pions, the invariant mass of at least one kaon pair is required to 
be within ±12 MeV/c 2 of the nominal value of the meson mass Ej. This requirement 
means that the D meson sample is dominated by decays containing an intermediate 
meson. A fit requiring the final state particles to originate from a common point is made 
and the \ 2 per degree of freedom (x 2 /ndf) of this fit is required to be less than five. In 
order to remove poorly reconstructed candidates, a cut is made on the uncertainty of the 
reconstructed invariant mass estimated by propagation of the individual track covariance 
matrices. The value of this cut depends on the decay mode under consideration and 
is chosen such that the bulk of the distribution is kept and only events in the tail are 
rejected. In a few percent of the events the reconstruction procedure gives rise to duplicate 
candidates. Therefore, if two or more candidates that are separated by less than 0.05 in 
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pseudorapidity and 50 mrad in azimuthal angle are found within one event, only that with 
the best D vertex x 2 is kept. 

Each candidate D meson, selected in this way, is combined with a well-identified muon 
that is displaced from the pp interaction vertex (impact parameter x 2 > 4) and that has 
Pt larger than 800 MeV/c to form a B candidate. A fit is made requiring the muon and 
the D candidate to originate from a common point and the x 2 P er degree of freedom of 
this fit is required to be less than five. To select semileptonic B decays, the invariant mass 
of the B candidate is required to be in the range 2.5 — 6.0 GeV/c 2 . In principle, the large 
combinatorial background can be further reduced by cutting on the decay time of the 
D meson, but due to the correlation between the decay time and the mass, this cut would 
bias the mass distribution. Therefore, a cut requiring significant displacement between the 
6-hadron decay vertex and the associated pp interaction vertex is applied. This achieves 
high signal purity whilst not biasing the distribution of the D decay time. 



4 Fit Results 

The D meson masses are determined by performing extended unbinned maximum likelihood 
fits to the invariant mass distributions. In these fits the background is modelled by an 



exponential function and the signal by the sum of a Crystal Ball 33 and a Gaussian 
function. The Crystal Ball component accounts for the presence of the QED radiative tail. 
Alternative models for both the signal and background components are considered as part 
of the studies of the systematic uncertainties. The model for the signal shape contains six 
parameters: 

• a common mean value for the Gaussian and Crystal Ball components; 

• the widths of the Gaussian (ctq) and the Crystal Ball (ctcb) components; 

• the transition point (a) and exponent (n) of the Crystal Ball component; 

• the relative fraction of the Crystal Ball (/cb) component. 

To reduce the number of free parameters in the fit, a, n and /cb together with the ratio 
of o"cb to ctq, are fixed using a simulation that has been tuned to reproduce the mass 
resolution observed in data for the B + — > J/ipK + and B + — > J/ipK + 7r~7r + decay modes. 
By fixing the ratio of o"cb to a G the resolution model is constrained up to an overall 
resolution scale factor that is close to unity. 

The Crystal Ball function describes the effect of the radiative tail far from the peak 
well. However, close to the peak its shape is still Gaussian, which results in a bias on the 
fitted mass that scales with the Q-value of the decay mode. This effect is studied using 



Photos 27 to model the effect of QED radiative corrections. The size of the bias is found 
to be 0.03 ± 0.01 MeV/c 2 for the D° -»■ K + K-K~tt + mode. For the D° -»■ K+R-n+n-, 
D + K + K~n + and D+ K + K~ir + decay modes a value of 0.06 ± 0.01 MeV/c 2 is 
found. These values are used to correct the mass measurements. The effect cancels in the 
measurement of the mass differences. 
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Figure 1: Invariant mass distributions for the (a) K + K~ir + ir~ and (b) K + K~K~ir + final states. 
In each case the result of the fit described in the text is superimposed (solid line) together with 
the background component (dotted line). The pull, i.e. the difference between the fitted and 
measured value divided by the uncertainty on the measured value, is shown below each plot. 

Table 1: Signal yields, mass values, resolution scale factors and binned x 2 / n df (using 100 bins) 
obtained from the fits shown in Fig. [T] and Fig. [2] together with the values corrected for the effect 
of QED radiative corrections as described in the text. 



Decay mode 


Yield 


Fitted mass 
[MeV/c 2 ] 


Corrected mass 
[MeV/c 2 ] 


Resolution 
scale factor 


X 2 /ri d f 


L> u — > K+K-n+n- 
D° -> K+K-R-TT+ 
D+ -> K+R--IT+ 
D+ -> K+K-JT+ 


4608 ± 89 
849 ± 36 
68, 787 ±321 
248, 694 ± 540 


1864.68 ±0.12 
1864.73 ±0.15 
1869.44 ± 0.03 
1968.13 ±0.03 


1864.74 ±0.12 

1864.75 ±0.15 
1869.50 ±0.03 
1968.19 ±0.03 


1.031 ±0.021 
0.981 ±0.042 
0.972 ±0.003 
0.971 ±0.002 


0.83 
0.92 

2.5 



The resulting fits for the D° decay modes are shown in Fig. [j] and that for the K + K~n + 
final state in Fig. [2} The values obtained in these fits are summarized in Table [1} The 
resulting values of the D + and masses are in agreement with the current world averages. 
These modes have relatively large Q-values and consequently the systematic uncertainty 
due to the knowledge of the momentum scale is at the level of 0.3 MeV/c 2 . Hence, it is 
chosen not to quote these values as measurements. Similarly, the systematic uncertainty 
due to the momentum scale for the D° — » K + K~tt + tt~ mode is estimated to be 0.2 MeV/c 2 
and the measured mass in this mode is not used in the D° mass determination. 

The quality of the fits is judged from the x 2 /ndf, quoted in Table[TJ and the fit residuals. 
It has been checked using simulated pseudo-experiments that the sizeable trends seen 
in the residuals for the K + K~tt + mode, where the dataset is largest, do not bias the 
mass difference measurement. The fitted resolution scale factors are all within a few 
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Figure 2: Invariant mass distribution for the K + K~n + final state. The result of the fit described 
in the text is superimposed (solid line) together with the background component (dotted line). 
The pull, i.e. the difference between the fitted value and the measured value divided by the 
uncertainty, is shown below the plot. 



percent of unity, indicating that the calibration parameters obtained from the B + study 
are applicable in this analysis. The uncertainties on the masses reported by the fits are in 
good agreement with the results obtained in pseudo-experiments. 

Using the values in Table [TJ the mass differences are evaluated to be 

M(D+) - M(D°) = 4.76 ±0.12 (stat) MeV/c 2 , 
M(D+) - M(D+) = 98.68 ±0.03 (stat) MeV/c 2 , 

where the uncertainties are statistical only. 



5 Systematic uncertainties 

To evaluate the systematic uncertainty, the complete analysis is repeated, including 
the track fit and the momentum scale calibration when needed, varying within their 
uncertainties the parameters to which the mass determination is sensitive. The observed 
changes in the central values of the fitted masses relative to the nominal results are assigned 
as systematic uncertainties. 

The dominant source of uncertainty is the limited knowledge of the momentum scale. 
The mass fits are repeated with the momentum scale varied by ±0.03%. A further 
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Table 2: Systematic uncertainties (in MeV/c 2 ) on the mass measurements and on their differences. 



Source of uncertainty 


M{D°) 


M(D+) - M(D°) 


M (£>+) - M{D + ) 


Momentum scale 


0.09 


0.04 


0.04 


Energy loss correction 


0.03 


0.06 


<0.01 


mass 


0.05 


<0.01 


<0.01 


Signal model 


0.02 


<0.01 


<0.01 


Background model 


<0.01 


<0.01 


<0.01 


Quadratic sum 


0.11 


0.07 


0.04 



uncertainty is related to the understanding of the energy loss in the material of the 
tracking system. The amount of material traversed in the tracking system by a particle is 
known to 10% accuracy [34]. Therefore, the magnitude of the energy loss correction in 
the reconstruction is varied by ±10 %. 

Other uncertainties arise from the fit model. To evaluate the impact of the signal model, 
a fit is performed where all signal parameters are fixed according to the values found in 
the simulation and a second fit where the parameters ctq and <tcb axe allowed to vary 
while keeping the relative fraction, /cb, of the two components fixed. The larger of the 
differences to the default fit result is assigned as an estimate of the systematic uncertainty. 
Similarly the effect of the background modelling is estimated by replacing the exponential 
function with a first-order Chebychev polynomial. The shifts of the mass values observed 
in these tests are generally much smaller than 0.01 MeV/c 2 and are assigned as systematic 
uncertainties. For the K + K~tt + fit further cancellation occurs in the mass difference. It is 
concluded that the details of the fit model have little effect on the presented measurements. 

An additional uncertainty arises from the knowledge of the value of the K + mass, 
rriK± = 493.677 ± 0.016 MeV/c 2 (4). The effect of this uncertainty on the measurements 
has been evaluated from simulation studies. 

The systematic uncertainties on the measured masses and mass differences are sum- 
marized in Table [2] The uncertainties related to the momentum scale and energy loss 
correction are fully correlated between the measurements. 

Various cross-checks of the measurements are made. Two checks are related to the 
knowledge of the tracking system alignment. First, a study has been performed where 
particle trajectories are reconstructed without using the information related to the tracking 
detector located before the entrance of the spectrometer magnet. This information is not 
required to form a track but improves the momentum resolution by 10 — 20 %. The second 
test is to vary the track slopes in the vertex detector by the uncertainty of 2 x 10~ 4 on the 
length scale of the detector described in Ref. 35 . The results obtained in these studies 



are consistent with those presented here and no additional uncertainty is assigned. 

A further check for the D° mass measurement is the comparison of the measured mass 
in the D° — > K + K~K~ir + mode with that obtained in the three other four-body modes. 
Systematic effects related to the momentum scale will affect modes with a high Q-value 
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more than those with low Q-values. The relationship between the reconstructed mass (to) 
and the momentum scale (a) after a first-order Taylor expansion in m 2 /p 2 is 



m 



m 



true 



a] 



where 



TO7 



Pi 



p is the total momentum of the decaying meson and and rrii are the momenta and 
masses of the daughter particles. This formalism assumes that there are no additional 
differences affecting the momentum scale between the modes such as differences in track 
kinematics or the effect of QED radiative corrections. For each decay mode the average 
value of / is obtained from the data using the sPlot technique [36] with the mass as the 
control variable to subtract the effect of background. The values obtained in this way are 
in good agreement with those found in the simulation. In Fig. [3] the measured D° mass is 
plotted versus / for the four-body decay modes studied here. The shaded area on this 
plot corresponds to the assigned systematic uncertainty of 0.03% on the momentum scale. 
Though there is evidence of a systematic effect for the low /-value modes it is accounted 
for by the assigned uncertainty. 
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Figure 3: Measured D° mass versus / as defined in Eq. [5j The (yellow) shaded area corresponds 
to a systematic uncertainty on the momentum scale of 0.03 % centred on the result for the 
D° — > K + K~ K~tt + mode (horizontal dashed line). Only the D° — > K + K~ K~ir + mode, where 
the systematic uncertainty is lowest, is used to determine the D° mass. 
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Table 3: LHCb measurements, compared to the best previous measurements and to the results 
of a global fit to available open charm mass data. The quoted uncertainties are the quadratic 
sums of the statistical and systematic contributions. All values are in MeV/c 2 . 



Quantity 


LHCb 


Best previous 






pdg fit g 


measurement 


measurement 








M(D°) 


1864.75 ±0.19 


1864.85 ±0.18 


5 




1864.86 ±0.13 


M(D + ) - M(D°) 


4.76 ±0.14 


4.7 ±0.3 


7 




4.76 ±0.10 


M(D+) - M(D + ) 


98.68 ±0.05 


98.4 ±0.3 


10 




98.88 ±0.25 



The dataset has also been divided according to the magnet polarity and data-taking 
period and for the charged modes by the sign of the product of the magnet polarity and 
the D meson charge. In addition, for modes where the event samples are sizable the 
measurements are repeated in bins of the D meson kinematic variables. None of these 
tests reveal any evidence of a systematic bias. 

6 Summary 

Measurements of D meson masses and mass differences have been performed using pp 
collision data, corresponding to an integrated luminosity of 1.0 fb _1 collected at a centre- 
of-mass energy of y's = 7 TeV with the LHCb detector. The results are 

M(D°) = 1864.75 ±0.15 (stat) ±0.11 (syst) MeV/c 2 , 

M(D + ) - M(D°) = 4.76 ±0.12 (stat) ±0.07 (syst) MeV/c 2 , 
M(D+) - M{D + ) = 98.68 ±0.03 (stat) ± 0.04 (syst) MeV/c 2 . 

The dominant systematic uncertainty is related to the knowledge of the momentum scale. 

As shown in Table [3j these measurements are in agreement with previous measurements. 
The results for the mass differences have smaller uncertainty than any previously reported 
value. The measured value of the D° mass has a similar precision to the published CLEO 
result [5]. Including this result in the determination of the X(3872) binding energy given 
in Section [l] gives Eb = 0.09 ± 0.28 MeV/c 2 . This reinforces the conclusion that if the 
X(3872) state is a molecule it is extremely loosely bound. 

The measurements presented here, together with those given in Ref. |4| for the D + and 
D° mass, and the mass differences M(D + ) - M(D°), M(D+) - M(D + ) can be used to 
determine a more precise value of the mass 

M(D+) = 1968.19 ± 0.20 ± 0.14 ± 0.08 MeV/c 2 , 

where the first uncertainty is the quadratic sum of the statistical and uncorrelated system- 
atic uncertainty, the second is due to the momentum scale and the third due to the energy 
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loss. This value is consistent with, but more precise than, that obtained from the fit to 
open charm mass data, M(D+) = 1968.49 ± 0.32 MeV/c 2 g). 
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